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Abstract 

The effect of correlated scalar-isoscalar two-pion exchange (CrTPE) modes is considered in con- 
nection with central and spin-orbit parts of the AW force. The two-pion correlation function is 
coupled directly to the scalar form factor of the nucleon which we calculate in the large- N c limit 
where the nucleon can be described as a soliton of an effective chiral theory. The results for the 
central AW force show a strong repulsive core at short internucleon distances supplemented by a 
moderate attraction beyond 1 fm. The long-range tail of the central AW potential is driven by 
the pion-nucleon sigma term and consistent with the effective a meson exchange. The spin-orbit 
part of the AW potential is repulsive. The largeW c scaling behavior of the scalar-isoscalar AW 
interaction is addressed. We show that the spin-orbit part is 0(1/N^) in strength relative to the 
central force resulting in the ratio ~ 1/9 suggested by the 1/N C expansion for N c = 3. The latter 
is in agreement with our numerical analysis and with the Kaplan-Manohar large- A,; power count- 
ing. Unitarization of the tttt scattering amplitude plays here an important role and improves the 
tree level results. Analytical representations of the CrTPE AW potential in terms of elementary 
functions are derived and their chiral content is discussed. 
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I. INTRODUCTION 



Understanding the intermediate and short range parts of the nucleon-nucleon (NN) interaction 
is still an interesting problem and much effort is being put into this topic both from theoretical 
and from experimental sides. The long-range part of the NN interaction is well described and 
represented by means of the one-pion-exchange (OPE) potential. For higher momentum transfer 
or small internucleon distances the NN dynamics becomes complex and one has to rely on phe- 
nomenological models. In this region the a meson [lj], for which the empirical evidence remains 
controversial, effectively represents scalar-isoscalar multi-pion correlations and generates an inter- 

and one-boson-exchange 



iiion 

3,3 



mediate range attraction in microscopic QCD motivated approaches 

(OBE) models ], |J. The latter are most popular and phenomenologically successful in describing 
the nuclear force. 

With the recent progress of chiral perturbation theory (%PT) in meson-meson and meson-baryon 
sectors the idea to extend it to the description of the NN interaction appears to be quite natural. 
The chiral SUl(Nf) x SUn(Nf) symmetry provides a consistent framework for the construction 
of the NN force HQ, [J. In ^PT the NN potential has already been calculated in N 3 LO H 1(| 



and much wor 
information 



c has been done in applying these ideas to nuclei, too. Supported by empirical 



111 ], the calculations are important and transparent for the peripheral NN partial 



11 



waves, which are well understood in terms of OPE and two-pion-exchange (TPE) processes 
ll^ |. The long-range behavior of the TPE force and the effect of correlated cr-like scalar-isoscalar 
7T7T modes are related problems. The su ccess a nd dynamical content of OBE models suggest that 
meson-meson correlations are important 



ill 

This fact is motivated by dispersion theory, where 



the main effects of higher-order interactions can be accounted for by the inclusion of experimentally 
known meson resonances, which bring singularities of the amplitude close to the physical region. 
In the dispersion-theoretical framework, the a can be explained by correlated two-pion exchange 
with a broad spectral distribution around ~ 4:M^ which ultimately leads to the isoscalar central 
attraction between two nucleons. 

In baryon xPT the tttt correlations are taken into account at tree level and the relation to the 
genuine a exchange is lost. For example, in heavy baryon xPT the exchange of two uncorrelated 
pions, formulated using chiral symmetry constraints and including A isobars, explains the tail 
of the scalar-isoscalar NN potential without any need for a true scalar meson |l5(. Additional 



consideration of tree level tttt correlations leads to the surprising result, that these terms are 
very small and - even more - lead to a weak repulsion Note, that in the effective theory 

without explicit A degrees of freedom the chiral invariant NN interaction terms are accompanied 



in subleading order by low-energy constants 
central potential is dominated by c\ and C3 



LEC) ci,C3 and C4 |l6(]. In particular, the isoscalar 



111 ] proportional to the pion-nucleon sigma term 



and nucleon axial-polarizability, respectively. The LEC C3 can be approximately saturated by 
the A resonance. But the dimension-two operators q include also information about i-channel 
meson exchange suggesting a strong influence of i-channel singularities. The latter was studied in 
Refs. 17, 15 ] (see also Ref. where it was shown that the exchange of a scalar meson with mass 
and coupling constant similar to a allows to explain c\ and, interestingly, also part of c%. Finally, 



C4 mentioned above is dominated by isovector tttt correlations from p meson exchange [l7j. Since 
we concentrate in this paper on s-wave CrTPE, i.e. a- meson channel, in the following we will 
discuss matters predominantly connected with c\. 

As well known the tree level contact interaction does not account for the entire tttt dynamics, and 
possibly one has to go beyond the tree level approximation to restore the relation to t 
Recently, non-perturbative methods describing tttt scattering were developed 



ne a meson. 



20, 



21 



22, 



23, 



24]. 



One of them was proposed in Ref. 25], where based on the Bethe-Salpeter description of the 
tttt scattering process the re-summation of the infinite series of pion loops was suggested. This 
method which only uses the lowest order chiral Lagrangian as input and implements an exact 
unitarity, leads to a dynamical pole in the tttt amplitude with a position around M — iT/2 ~ 
450 — i221 MeV. Naturally, it is identified with the a meson. Considering the role played by the 
unitary tttt correlations (or dynamical a state) in the NN interaction problem, it was found [3] 
that the resulting scalar-isoscalar central NN potential shows a strong repulsion at internucleon 
distances less then ~ 1 fm, and a moderate attraction at r > 1 fm. This feature of unitary 
tttt scalar-isoscalar correlations differs from the conventional Yukawa- like a exchange, V a (r) ~ 
— exp(— m a r)/r, which always results in an attraction. Interestingly, in phenomenological NN 



2J and CD-Bonn 



301 the isospin 



3J|. But the 



interaction models like Paris 27], Argonne V18 [28j |. Nijmegen 

independent scalar components are rather similar for AW separations beyond 0.5 fm 
radial dependencies differ considerably at short distances, ranging from attractive for Nijmegen and 
CD-Bonn (by construction) to repulsive for the more phenomenological V18 and Paris potentials. It 
also was argued that the appearance of a strong repulsion in the scalar-isoscalar channel cannot be 



interpreted in terms of meson exchan 



a remarkable feature which only was known from the 



Skyrrno model |^ and now from R<\T. 2(>'. [t is instructive. I hat in the Skvnno model the use of lire 
simplest product ansatz results in a repulsive central force [33( , and the missing intermediate range 



34] or additionally 



35]. 



attraction can be produced only if one goes beyond the product approximation 
considers the correlated two-pion exchange (CrTPE) in the scalar-isoscalar channel 

In this work the role played by the CrTPE in the NN interaction is reconsidered for the con- 
struction of the scalar-isoscalar central and spin-orbit NN potentials. Considering the tree level 
and unitary two-pion correlations in the NN interaction problem, we couple the 7T7T correlation 
function to the scalar form factor of the nucleon in a model-independent way. The latter is cal- 
culated in the limit of a large number of colors (large-iV c limit) using the Chiral Quark-Soliton 



Model (xQSM) 



36 



where the nucleon can be described as a soliton of an effective chiral theory. 



This way our results support the structure of the scalar-isoscalar central NN force observed in the 



unitarized xPT 



and in the Skyrme model 



341 ] . We find a strong repulsion at short distances 



and a moderate attraction at r > 1 fm - even with the tree level 7T7T interaction. The tail of the 
central CrTPE potential is driven by the the pion-nucleon (ttN) sigma term and is consistent with 
the effective a exchange in OBE models. By this, the effective oNN coupling constant can be 
related to the ttN sigma term. We also show that the above feature of the central NN force is 
general. It mainly relies on the particular functional form of the scalar form factor of the nucleon 
and is not the effect of unitarization. Another point of interest in this work is the generation of 
the scalar-isoscalar spin-orbit (LS) force. Our results for this part of the NN potential show a 
repulsion - a feature which differs from phenomenological a exchange resulting in an attractive LS 
interaction. 

The large- iV c behavior of CrTPE forces is considered and consistency with large- N c QCD anal- 
yses and the large- N c behavior of scalar-isoscalar components of phenomenological NN interaction 
models is shown. Both, the central and spin-orbit parts satisfy large- N c QCD counting rules. We 
find that the relative strength of spin-orbit and central interactions scales according to ~ 1/9 in 



3J. Addi- 



remarkable agreement with the Kaplan-Manohar large- N c spin-flavor power counting 
tionally, analytical representations of central and LS potentials in terms of elementary functions 
are derived and their chiral content is discussed. We show that the chiral symmetry breaking part 
of the interaction is small and that terms non-vanishing in the chiral limit are dominant and fully 
drive all essential features of the CrTPE force in the scalar-isoscalar channel. 



The paper is organized as follows. In the next two sections we discuss the construction of the 
CrTPE scalar-isoscalar NN force. We emphasize the importance of the scalar form factor of the 
nucleon and give its phenomenological description in Sec. V. The large- iV c scaling and consistency 
are addressed in Sec. IV. Finally, the analytical results and discussions are presented in Sec. VI. 

II. TWO-PION CORRELATION FUNCTION 

We start our considerations of the scalar-isoscalar CrTPE NN interaction from the standard 
definition of the tree level tttt correlation function and its large-Ac behavior. To lowest order in 

(2) 

the derivative expansion the mesonic Lagrangian CX is given by the SU (3) nonlinear a model 

f2 f2 r , 

£ (2) = ZjLTr^/^t^f/j +B ^L Tr ^([/ + [ /t)j (!) 

and contains the most general low-energy interactions of the pseudo-scalar meson octet. In Eq. Q 
the leading symmetry-breaking term is linear in the quark masses A4 and is characterized by 
B = —(qq)/f%. The pion-decay constant /„- = 93 MeV and the scalar condensate (qq) are the 
relevant parameters. In the SU(2) limit the non-linear field U entering Eq. Q is given by the 
standard matrix form: 

U(x) = exp (»*(*)//*) , $ = T7T (2) 

where r are the SU(2) Pauli matrices and it is the isovector pion field. For the ir a irb — ► n^d 
scattering process, defined by the Cartesian isospin indices a,..., the use of the standard xPT 
procedure in expanding the jfifr to order 0(tt a ) results in the tree level contact interaction 

- iV£T cd = S ab 5 cd A(s) + 5 ac 5 bd A(t) + S ad 5 bc A(u), (3) 

where 

^^^) = j2^~Ml- \ E A l )+0(q% (4) 

71 i=a,b,c,d 

and Aj = kf - Ml are the off-shell part of the invariant tttt amplitude. The Mandelstam variables 
are related by s + t + u = k\ + + k\ + k\. At this order of the pion field expansion the isoscalar 
5-wave tttt partial amplitude (L = 0) is obtained from the standard decomposition 



where Pl(cos9) are the Legendre polynomials and (V2) a accounts for the statistical factor oc- 
curring in states with identical particles: a = 2 for tttt — > tttt. The tree level scalar-isoscalar tttt 



scattering amplitude is 



V, 



L=I=0 



M 2 
~2 



(6) 



f 2 V 2 3 

In the limit of a large number of colors N c of QCD, the scaling behavior of Eq. © can be 
obtained from large- N c xPT 39]. The conventional large- iV c counting rules require that the pion 
mass is M n ~ 0{1) and the decay constant is f v ~ 0{y/W c ). As a result the tttt amplitude scales 
according to 

V m ~ 0(l/N c ) (7) 
and large- N c QCD becomes a theory of weakly interacting pions. Note that the unitarization of the 



■kit correlation function (scattering amplitude) V 7 



l=i=o 



Eq. ©, is a 



23j and Bethe-Salpeter 



24 



engthy procedure. Here, we 



3, 



only mention some recent analyses which employ dispersion relations 21], Pade approximants 
solution of Roy 



22], 



provided by most of these methods contain a dynamical pole 



25] equations. Unitary tttt scattering amplitudes as 
{20! which is identified with the a 



meson. 



III. TWO-PION EXCHANGE KERNEL AND SCALAR FORM FACTOR 



With Eq. © we come to the definition of the CrTPE potential in the scalar-isoscalar channel 
describing the on-shell NN scattering process, N(pi, si)N(p 2 , S2) — ► ^(^3,53)^(^4,54), where Pi 
and Si are the four momenta and spins of interacting nucleons, respectively. As well known, in 
the generic case of TPE the chiral symmetry constraints work best and allow to relate the TPE 
process containing two intermediate pions, with Cartesian isospin indices a and b to the off-mass 



shell ttN scatter ing amplitude in a model-independent way 



Green function 



8, 



10 



3 



. The latter is defined by the 



N 



d 4 x d 4 y (p',s'\fP a (x)P b (y)\p,s)e i{k ' x - k ^ 



Tf N (p',k',p,k) 



(8) 



k 2 -M% k' 2 -M 2 

Here G w accounts for the coupling of the pion field to the pseudoscalar quark densities 
(0|P a (0)|7r b ) = 5 a bG n , P a (x) = i^(x)'j5T a ilj(x) and p',p (k',k) are the momenta of outgoing and 
incoming nucleons (pions). For the irN scattering amplitude the standard decomposition 
involves four Lorentz invariant functions D^^B^ : 



J jab 
■kN 



5ab F+ N + ^[T a ,T b ] F nN , 



(9) 



o 



FIG. 1: The coupling of the scalar-isoscalar tttt correlation function V T 



l=i=o 



symmetric part F^ N (filled disk) of the ttN scattering amplitude. 



F ± 



u(pj 



AT) ' 



D ± + ia» u k' ll k u B ± 



u(p N ,s). 



(filled box) to the isospin 



(10) 



We are interested in the chiral content of the TPE potential in the scalar-isoscalar channel where 
the introduction of the tree level contact 7T7T interaction results in a NN scattering amplitude with 
two connected pion loops 



IVnn 



d A k d A k 



(2tt) 4 (2tt) 4 



V, 



£=7=0 



(11) 



(k 2 - M 2 ){k> 2 - M 2 ) 



(k 2 - M 2 ){k' 2 - M 2 ) 

where k{k) and k'(k') = k{k)+p' —p are the momenta of the exchanged pions, p{p') are relative NN 
momenta in the initial (final) state and the superscripts refer to nucleons (1) and (2), respectively. 
F^ N is the isospin symmetric part of the ttN amplitude, Eq. Q, and the scalar-isoscalar tree level 
contact interaction V^ I=0 is given by the Eq. Q. The diagrammatic representation of Eq. 
is shown in Fig. ^ 

The above expression, Eq. 1)1 ljl has well known difficulties. It involves the off-mass-shell tttt 
interaction which is not unique and depends on the choice of the interpolating pion field (parame- 



terization of the SU (2) matrix U) 



15, 



431 ] . If the initial and final pions are all on the mass shell, 
i.e. Aj = in Eq. ©, the tttt scattering amplitude is independent of a change of field variables in 
agreement with the equivalence theorem ^ | . In other words Eq. © gives a unique result indepen- 
dent of the parameterization of U only for the on-shell matrix elements. To avoid these ambiguities 
it was stated in Ref. 15] that one has to include a subset of diagrams to find cancellations of the 
off-shell tttt isoscalar amplitude. This statement was rigorously verified in Ref. j^fj], where it was 
shown that the additional consideration of a subset of chiral diagrams results in exact cancellations 
of the off-shell part of the tttt amplitude and the on-shell residue can be factorized out from the 
loop integrals of Eq. (|llj) . By this and using notations of Ref. [2] the NN amplitude in the general 



case of L = I = tttt exchange takes the form 

v NN (t) = 6 [r+(t)]( 1 )[r+(t)]( 2 ) v^=°(f), (12) 

where V^ r =/=0 is the on-mass shell part of Eq. (JJJJ) with s — ► i = (p 7 — p) 2 = — g 2 in the A~A" cm. 
frame. In Eq. (|12[) the "vertex functions" [Ti-(t)] are defined by 

H ^ l jJ " 2 7 (2tt) 4 (fc 2 - M^)(F - M 2 ) 1 ^ 

and can be interpreted in the heavy and relativistic baryon xPT as one-pion-loop contributions to 
the scalar form factor of the nucleon a(t) [lH 4fl |: 

[r^W] (l) = ^[H (i) - (14) 

At t = the (j(0) is referred to as the pion-nucleon sigma-term, a n N. Note, that by construction 
the factor ~ l/-^ 2 entering Eq. (|14|) is a direct reflection of the chiral symmetry breaking part 
of the xPT Lagrangian, Eq. (|T|). The resulting AW scattering amplitude can be summarized as 
follows 

v NN (t) = I (^|) 2 M (1) M {2) v^ =/=0 (t) (15) 

The quasi-potential reduction of the AW amplitude Eq. ()15|) is a standard procedure. The 

Partovi-Lomon method is the most popular one and allows to reduce the Bethe-Salpeter kernel to 

I I 

the quasi-potential containing minimal relativity [46J]. The latter is ready for the iteration in the 



261 ] and wish to discuss the non-relativistic 



Lippmann-Schwinger equation. Here we follow Ref. 
potential only. With the normalization for Dirac spinors, u(p, s)u(p, s) = 1, the non-relativistic 
CrTPE potential is obtained from Eq. ()15j) simply by keeping its energy independent parts: 



2 



In square brackets the first term gives the central CrTPE potential Vc(t), with the structure similar 
to that derived in Ref. and the second term provides the additional spin-orbit part VLs(i) 
with the spin content: $7^5 = i(cr^ + cr^)(p' x p)/2. In coordinate space the central potential is 
given by the Fourier transform of Vc(t) and the spin-orbit (LS) part of the A" A" force is related to 



the central potential simply by 



46] 



1 l dV c (r) 

f I VLs{r) = -2M*r-dT' (17) 

where the spin-orbit operator £Ils m r-space representation has been omitted. The central part 
of Eq. (fTB|) has an interesting feature. Its long-range tail or the value of the potential at t = is 



determined by the pion-nucleon sigma term: Vc(0) = a^ N /3M£f£. So, the on 



unknown element 



12,45] Eq. (II) and 



which enters Eq. (|16f) is the scalar form factor of the nucleon. As noted in Refs. 
therefore also Eqs. (|15|) and are general and should be independent of models or approximation 
schemes used to calculate the scalar form factor of the nucleon. In this work, we calculate a(t) in 
the large- N c limit in the framework of the xQSM. But before doing model calculations the stated 
large-iV c behavior of Eq. 1)16(1 must be understood and consistency with large--/V c QCD analyses 
must be shown. 



IV. LARGE- iV c SCALING AND CONSISTENCY 

In this section we explore qualitative features of the scalar-isoscalar components of the CrTPE 
force Eq. (|16|) which may be understood directly from QCD. In this context, it may be useful to 



consider the interaction in the limit, when the number of colors, iV c , of QCD becomes large 



48], 



and to treat 1/N C as an expansion parameter. Some features of the NN force can be obtained in 



this limit in a model-independent way 



38, 



49 



501 ] . Recently, it was realized that large- iV c nuclear 



interactions are spin-flavor SU (4) symmetric, with Wigner's supermultiplet symmetry following as 
an accidental symmetry 49]. If we assume, in addition, that the spin-flavor symmetry properties 
of the NN interaction are independent of phases of the many body ground state (at N c — > oo 
nuclear matter forms a classical crystal, and a phase transition between N c = 3 and N c = oo is 
expected 38]), then the large- N c scaling behavior of the NN force can be analyzed in a general way 



38 



1 



resulting in QCD expectations which can explain general features of the NN interaction 

The first treatment of the NN interaction in large- iV c QCD was done in Ref. [3], where il was 
argued that the dominant interaction between two baryons is of order ~ 0(N C ). This expectation is 
consistent with nuclear dynamics where the dominant interaction components are isoscalar-scalar 
and -vector interactions 

V a ~ 0(N C ), V w ~ 0{N C ). (18) 
In phenomenological models, the latter can be parameterized in terms of effective a and u ex- 



changes 



The entire analysis of the NN potential was done in Ref. |38| where it was found 



that large- A^ c QCD implies that the scalar-isoscalar central and spin-orbit forces scale as 

V c ~0(iV c ), V LS ~0(l/N c ). (19) 



Recently, the dynamical interpretation of the NN interaction for four modern phenomenologically 



successful models 



28, 



29 



30] in the large- iV c limit was reported and consistency of their scalar- 



isoscalar components with Eq. (|19|) was found 

In our case, the consistency condition implies that the large-iV c scaling behavior of Eq. H16|) 
must be consistent with large- N c QCD analysis p^ - with the large- N c scaling of scalar-isoscalar 



D 



components of all known phenomenologically successful NN interaction models |3J| and also with 
effective a meson exchange, V a ~ 0(N C ). The dynamical quantities which enter the scalar-isoscalar 
CrTPE force, Eq. (|16|) . are the scalar form factor cr(t), the pion mass M n , the nucleon mass Mjy and 
the 7T7T scattering amplitude V W7r which scales as V nn ~ 0(1/N C ), Eq (JJJ). We follow standard iV c 
counting rules |39j, |5J| where a(t) ~ 0(N C ), M w ~ 0(1), Mjy ~ 0{N C ) and the nucleon momenta 
|p|,|p'| ~ C(l)- Counting powers of iV c the central part of the CrTPE potential, first term in 
Eq. (|TT))) . is Vc ~ 0(N C ). Its spin-orbit counterpart VLsif) additionally contains two inverse 
powers of nucleon mass ~ ~ 1/N^ resulting in Vls ~ 0(1/N C ). By this, the consistency 

with scaling relations from large- iV c QCD Eq. (|19|) is noted. Interestingly, Eq. (1191) implies that 
the spin-orbit force Vls is 0{1/N^) in strength relative to the central force Vc [38j 

V L s/Vc ~ 0(l/N 2 c ) (20) 



We refer to Eq. ()20|) as the Kaplan-Manohar scaling relation [38j. In the last section it will be 
shown that Eq. ()20|) is very accurate and in remarkable agreement with our results, signaling that, 
indeed, with the actual number of colors N c = 3 the relative strength of potentials scales like 
Vls/Vc ^ 1/9. 



V. SCALAR FORM FACTOR IN THE LARGE-iV c LIMIT AND QUARK- 
SOLITON MODEL 



In the large-iV c limit, the QCD is equivalent to an effective theory of mesons with baryons 
emerging as solitonic configurations |48j. The xQSM 36| provides a practical realization of the 



large- A^c picture of the nucleon and is considered as a chiral relativistic quantum field theory of 
quarks, anti-quarks and Goldstone bosons. It is defined by the partition function in Euclidian 



space which is the path-integral over pseudoscalar meson and quark fields 52 
Z = j VipV^VU exp i J d 4 xi;(i Jd - MqlT 5 - m q )il> 



53] 



(21) 



where U denotes the SU(2) pion field, Eq. 



ET* = exp(i75*// ff ) = \{U + rf) + \{U - t/t) 75 



(22) 



In Eq. ()21() M g and m g are dynamical, arising from spontaneous breaking of chiral symmetry and 
current quark masses, respectively. The action, Eq. (|21|) . was derived from the instanton vacuum 
model |52j |. where the cut-off is determined by the average size of instantons (p), and the dynamical 
quark mass M q is momentum dependent. In the large- N c limit the chiral bosonic fields U" /5 can be 
integrated out by the saddle point method using a classical background field of hedgehog shape 



(23) 



which plays the role of a Hartree-type mean field for quarks forming a soliton-like bound state. 
With r = \x\ and r = x/r, the variational procedure reduces to the determination of the self- 
consistent soliton profile V(r), where V(r — ► oo) = and "P(O) = —tt produce a soliton with unit 
winding number. In practical calculations the soliton size R s = 1/M q is treated as a free parameter, 



playing the role of the axial coupling constant gA 



n 



lim 

r— >oo 



lim r V(r) 

m,r— >0 



-2Ri 



(24) 



Nucleon states of definite spin and isospin are obtained by quantizing the rotational zero modes of 
the soliton. 

The model expression for the scalar form factor of the nucleon in the xQSM is quadratically 
ultraviolet (UV) divergent and requires regularization 54 , 



55] 



a(t) = m q N c / d 3 xj { 



tx 



reg 



(25) 



Here T> n {x) = $* (x)"fQ$ n (x) is the scalar-quark density in the nucleon and jo is the spherical 
Bessel function. The completeness of quark Fock states | n) in the external pion field require that 
V n is represented as a sum over occupied (valence) and non-occupied negative energy Dirac-sea 
states. So it is useful to consider the contributions from discrete levels (lev) and from the Dirac- 
sea (sea) continuum separately. o~i ev (t) is finite and requires the single-particle wave functions 
\$>l ev (x)) to be found from the Dirac equation in the external pion field. The latter can be solved 
numerically. o~ sea (t) is UV divergent and depends on the regularization scheme employed to make 
it finite. To obtain o~ sea (t) we follow the method developed in Ref. [37 1. In short, the procedure is 
as follows. 



11 



To evaluate cr sea (i), the model expression for the continuum contribution is expanded in a series 



of the [/-field gradients - the interpolation formula method 



divergent, and (1 V-Jinik* parts. As was shown in Kef. [37 1 the latter is strongly suppressed with 
respect to the UV-divergent term by the instanton packing fraction or parametric smallness of 
M2(p) 2 < 1, and can be neglected. The dynamical quark mass is momentum independent and its 
value M q = 350 MeV is taken at zero momentum transfer from instanton phenomenology |57j |. We 
recall that the xQSM is defined with some appropriate regularization. In Ref. the regularization 
procedure is solved in a model independent way by observing that the structure of divergencies in 
a sea are the same as in model expressions for the vacuum scalar-quark condensate (ipip) and pion 
decay constant f n . The latter are fixed by their empirical values. The "arctan-profile" function 
which provides an accurate representation of the self-consistent profile is used 



561 ] . The series in VC/ contains UV- 



V{r) 



-2 arctan 



5f (1 + M w r)e 



-M n r 



Above steps result in a model expression for <J sea (t): 

/•OO 

a sea (t) = Ml ft / d 3 rj (i 
Jo 



t) 



2e 



(26) 



(27) 



where £ = (R^/r 2 )(l + M n r)e~ MlTr . It follows that due to the partial cancellation of additional 
contributions from discrete levels and second order continuum terms in the VJ7 expansion, the 
accuracy of Eq. (|2*Tj) in the limiting case a sea {t) — > a(t) is 0(15%). The latter corresponds to 
the extreme case that the nucleon scalar quark density is formed exclusively by the Dirac-sea 
(pion cloud). With the expected accuracy Eq. (|2?|) reproduces all features of the scalar form 
factor observed in exact xQSM calculations and also in %PT and compares well to lattice QCD 
results. Eq. (|27|) results in pion-nucleon sigma term a n N = <t(0) ~ 68 MeV, which is consistent 



with empirical information 



58 



this method. We refer to Ref. 



59]. The analytical simplicity of Eq. (|27|) is a big advantage of 



371 ] for details and use Eq. (|27f) for our numerical calculations. In 



addition, consider the pion mass expansion of the profile function V(r), Eq. (|26|) 



P(r,M OT 



-2 arctan 



(28) 



In the xQSM the soft pion limit, V{r^M v — > 0), is known to approximate well both the self- 
consistent profile and Eq. (|26|) . Following the method of Ref. we find that in this limit a(t) 
takes an even simpler form 



a(t) = ^ MlflRl j (R s ^tj2) e- R ^. 



(29) 
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If one would expand Eq. (|27|) in Taylor series around M n = 0, then Eq. Q29[) corresponds to the 
leading analytic ~ M% contribution to the scalar form factor or to the ttN sigma term, a^jy, defined 
by Eq. (|27|) . The quality of this representation will be discussed later. We only note that, a simple 
form of Eq. (|29j) is very successful, at least for our considerations, because it allows us to express all 
our final results for the r-space CrTPE potential in terms of elementary functions. Furthermore, 
Eq. ()29|) accounts for all details of the CrTPE force observed here and compares well with direct 
numerical calculation obtained with Eq. (|27j) . 

VI. RESULTS AND DISCUSSION 

A. The central CrTPE force: tree-level vs unitary tttt scattering amplitude 

In the discussion of our results we start from the scalar-isoscalar central NN potential in 
coordinate space, which is given by the Fourier transform of the first term in Eq. ()16|) . Our results 
for Vc(r) with the scalar form factor a(t) defined by Eq. ()27j) and the tree level tttt correlation 
function, Eq. @, are shown in Fig. [2] (dashed curve). The central potential Vc(r) shows a short- 
range repulsive core with a maximum at the origin, Vc{f = 0) ~ 600 MeV (see insert of Fig. [2J), 
and is supplemented by a moderate attraction with a minimum, Vc — —20 MeV, at intermediate 
distances r > 1 fm. This behavior is different from the conventional picture provided by the 
effective a exchange in OBE models and also from results obtained in the heavy baryon )(PT [rsl ]. 
In Ref. [15I the effect of the CrTPE in the scalar-isoscalar channel is very weak and repulsive. 
The latter means, that the t-dependence of the vertex functions or functional form of the scalar 
form factor is very important for determining the precise strength and behavior of the potential 
in all ranges of distances. At the same time one can find a remarkable agreement between our 
results and the analysis of Ref. 



261 ] . where the unitary tttt amplitude Eq. (|cS0|) has been coupled 
to the nucleon cloud at the one pion loop. In Ref. [3| the divergent loop- integrals were regulated 
by phenomenological vertex form factors. Being sensitive to the particular choice of the cut-off 
mass, the resulting scalar-isoscalar central NN potential shows a strong repulsion at internucleon 
distances less than 1 fm, and a moderate attraction of — 10(— 15) MeV, at r > 1 fm. 



To make our results comparable with that from Ref. 



261 ] we unitarize the tttt scattering ampli- 
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FIG. 2: The central CrTPE potential in coordinate space obtained with tree level (dashed curve) and 
unitarized (solid curve) rnr scattering amplitude. The dot-dashed and dash-dash-dotted curves correspond 
to the soft pion limit with M q = 350 MeV and M q = 410 MeV, respectively. The insert shows the entire 
structure of the NN potential. 



tude, Fig. |21 It can be done by the following substitution 



L=I=0 



Mi 



2\ - 1 



-i 



where G n7r (s) is a dimensionally regularized two-pion loop function (see e.g. Ref. |2fil 

1 



G n - !T (s) 



(47T) 



, , Ml , (7+1 

-1 + ln— £+<rln 



(30) 



(31) 



a — 1 

Here fj, = 1.1 GeV is the regularization mass fitted to the tttt phase shifts. The value of a in 
Eq. (|31j) is given by a - 



'\ _ -^l. Eq. ()31|) is analytic in space-like and time-like regions and 
for s > 4M% it develops an imaginary part, since a — 1 < and In = In + iir. For 
s — > t < the log behaves smoothly. Because Eq. (|30|) contains a pole in the s-channel around 
M n7T — iT/2 ~ 450 — i221 MeV, we restore the genuine relation to the cr-meson exchange, which now 
enters in our formalism as dynamical resonance in the tttt system. The main effect of the unitary 
tttt scattering amplitude (cr meson) in the t-channel is that it makes the repulsion ~ 400 MeV 
and attraction ~ — 10 MeV softer (solid curve in Fig. [2J). However it does not change the general 
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FIG. 3: Unitarization of the tttt scattering amplitude. 

features of the NN force already obtained at tree level. Interestingly, the latter effect is opposite 
to the role played by the the unitarization in the s-channel, where it leads to an enhancement of 
the interaction strength between pions. 

Another point is the relation of the effective a meson exchange of OBE models to the unitary 
scalar-isoscalar CrTPE force constructed here. In Ref. |3lJ the volume integrals were used to 
extract the strength of scalar-isoscalar components of phenomenological NN models and to relate 
them to the effective a exchange. It was shown that if the mass of the a is taken to be M a = 
600 MeV then the value of the aNN coupling constant for these interactions would range between 
9aNN = 7.6 (Paris), 9.0 (Argonne V18), 9.8 (Nijmegen) and 11.2 (CD-Bonn). It follows from 
Ref. 26], that the use of the unitary tttt scattering amplitude Eq. (|30|) results in an effective aNN 
coupling constant which can be expressed as 



9aNN 



V6V(t = 0){M*-M*/2) /U. (32) 



Here, V(t) is a vertex function which accounts for the coupling of the tttt correlation function to the 
nucleon pion cloud, and M a is an empirical a meson mass. Note that the use of M a = 600 MeV in 
this case is not entirely correct, because this value does not represent the actual pole of Eq. (|30|) . 
With the quoted value of V(t) at t = 26]: V(0) ~ 0.1 x 10~ 2 MeV^ 1 and the a mass (pole) 
M a ~ 450 MeV the resulting coupling is g a NN — 5, which has the right order of magnitude. In our 



case, note that the vertex function V(t) as provided by Ref. |26| is equivalent to our rt-(t) = V(t) 
in Eq. (|12f) . By this, we rewrite Eq. (|3*2*)) and relate the aNN coupling constant to the ttN sigma 
term a n N 

9*nn ^ yf^iv {Ml - /(f n M%). (33) 

Eq. 1(2*71) results in a n N ~ 68 MeV and our value for the V(t = 0) is given by: V ~ 0.12 x 
10~ 2 MeV -1 . With M a = 450 MeV our result is g a NN — 6.1. An artificial increase of the a mass 
to M a = 600 MeV, used in analysis of Ref. j^jj, results in g a NN — H which is just in the range of 
values from phenomenological models j^ . 

u 



Using the equivalence between T~^{t) and V{t) we inde pen dently verify the theoretical con- 
sistency of Ref. j^J. Indeed, with V(0) obtained in Ref. 



o^N — 57 MeV and is within the range of empirical values 
masses used in Ref. 



58, 



261 the resulting ttN sigma term is 



591 ] . It means that values of cut-off 



26j are reasonable and consistent with our results and with empirical infor- 



mation. There is also consistency between the large- N c behavior of Eq. (I33|) and large- N c QCD 
expectations, where a scalar-meson NN coupling constant is ~ 0{^fW c ) 38j- Because the 1/N C ex- 
pansion of the ttN sigma term follows the large- N c behavior of the scalar form factor a^jy ~ O{N ), 
the pion decay constant f n ~ 0{\/yfW c ), M w ~ 0(1) and M a ~ 0(1) [60j we get 



g aNN ~ 0(y/N c ), 



(34) 

in agreement with large- N c counting rules 3^ |. 

The noncommutativity of large- N c and chiral limits is a well-known phenomenon jfU]. The 
important physics behind this is the role of the A(1232) isobar. In the large- N c limit nucleon and A 
are degenerate with mass difference Ma — Mn ~ 0(1/N C ). As a consequence the contribution of A 
states is implicitly included in solitonic configurations and their effect on scalar-isoscalar quantities 
is twice that of nucleon states (37, 6^|. Note that, therefore the leading non-analytic contribution 
to a n N, which is 27 g\M% / (G<i-7r f%) obtained by expanding Eq. (|27j). is exactly three times larger 
than the corresponding value in xPT including nucleons and pions only. Being common for generic 
classes of hedgehog models 62|, this result can be explained by the presence of the intermediate A 



states in chiral loops. If one goes beyond the strict large- N c limit the effect of finite N — A mass 



splitting must be taken into account, the magnitudes of the resulting corrections are known 



62]. 



B. Saturation of the LEC c\ 



Here we use constraints imposed by chiral symmetry to check the consistency of the present 



1J. The 



approach. Consider the LEC c\ from xPT which is related to the ttN sigma term (J^n 
phenomenological interpretation of its value is based on the strongly coupled scalar-isoscalar meson 
exchange which can completely saturate ci, if 

l? " 180 MeV. (35) 



y/9aNN 

It is interesting to note that the effective a meson in the Bonn potential A\ with M a = 550 MeV 
and g^ NN /{An) = 7.1 respects this condition resulting in M a /^/g a NN — 179 MeV. In our case, 
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FIG. 4: The central CrTPE potential in the soft pion limit (solid curve). The parts non-vanishing and 
vanishing in the chiral limit are shown by the dashed and dot-dashed curves, respectively. 

using the pole of Eq. PU|I. i.e. M a = 450 MeV, and g a NN — 6.1, Eq. (J2HJ), the saturation condition 
is 

— g ~ 182 MeV. (36) 

y/9oNN 

which is in agreement with the value demanded from the scalar meson resonance saturation of the 
LEC ci. 



C. Analytical formulae for the central force 

The consideration of the central part can be continued in the soft pion limit where a{t) is given 
by Eq. (|29j) . In this case the central CrTPE force in the coordinate space representation can be 
expressed in terms of elementary functions. First, decompose the tree level tttt correlation function 
defined by Eq. © into two terms, Vf= I=0 = V$ + vj$: 

W = ~t/fl W=Ml/2fl (37) 
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By this, the CrTPE force reads as 



VNN(t) = V$ N (t) + V% ) N (t) 



(38) 



In the soft pion limit, due to the factor l/M% entering Eq. ()16j) . the first term V^ N does not 
contain any pion mass dependence and does not vanish in the chiral limit (non- vanishing part). 
The second term vffi N which follows vSk and contains an additional M% generates the symmetry 
breaking part of the NN force. With Eqs. (|29[) and (|37|) the Fourier transform of Eq. I)38|) can be 
carried out explicitly and after lengthy calculations (here we list our final results only) the central 
CrTPE force is given in coordinate space by following expressions 



^ flRs 



(39) 



80 + 8£ + 3£ 2 -3£ 3 



~\ 2 



2 + a 4 + £ + 2a/2£ 4 + £-2a/2£ 



Vg } (r,e<4) = 



Ml /4vr 2 



V^(r,e>4) = ^(^)/^ 




Ml / 2 V 2 ^ 

arctan \ \l — arctan 

2 I 2 \ 4-^ 




, ,e \ 1 1 / 2 v 2 ? 

;>■).'< '•( an I W — I — - < n - a re I an 




(40) 



(41) 



where £ = r 2 /R 2 and the superscripts (1) and (2) refer to the non-vanishing part and the symmetry 
breaking part, respectively. Note that Eqs. I)40j) and (|41|) are defined for £ < 4 and £ > 4, 
respectively. 

At small separation scales Eqs. ()39|) and (|4Uj) can be expanded in Taylor series around r — > 0, 
and up to C(r 4 ) the leading terms are given by: 



V?(r) 



vg\r) 



r->0 



V27T 2 



flMlRl 



3d - 



/ 2 M^r 2 + 0(r 4 ). 



(42) 



(43) 
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The strength of the central force is defined by the sum of Eqs. (|42|) and (|43j) at r = 0: Vc(0) 

5\/27r 2 ( ^ , 1,. 2d21 .,> . 



Vc(r) r=0 = l 1 + 5 r 8 > f * Rs - (44) 



For asymptotically large r —> oo (up to 0(r 8 )) the interaction is dominated by molecular-like 
NN forces 

V { c\r) = -64V27r 2 fiRl± + 0(r- 8 ), (45) 



vg\r) = *^MtfiRl\ (46) 

As seen from above expressions, and Vq do not vanish at small separations and both are 
repulsive. At asymptotically large r Vq becomes attractive, signaling that it should cross the 

(2) 

zero. Contrary, the symmetry breaking part, V c , stays repulsive for all distances and as a result 
of Vq /Vq ~ M^Rg/5 <C 1 it is suppressed relative to V^. The above qualitative arguments 
are supported by full analytical representations given by Eqs. (|39W41j) . The resulting central NN 
potential is shown in Fig. 0J It is clear from Fig. 0] that the part Vq \r), which is non- vanishing 
in the chiral limit (dashed curve), is dominant and fully drives the generic behavior of the central 

(2) 

CrTPE force (solid line). The chiral part V c (r), which is vanishing in the chiral limit (dot-dashed 
curve), is repulsive for all distances with a strength of about ~ 20 MeV at r = 0. Comparison with 
numerical results from the previous section reveals that at intermediate distances (dip region) the 
total CrTPE force obtained in the soft pion limit (shown by the dot-dashed curve in Fig. |5J) lies 
in-between the tree level and unitarized curves. At small separation scales (see insert of Fig. |5J) the 
integral Eq. (|27|) and analytical Eq. (|29|) representations for the scalar form factor work equally 
well, and the analytical and numerical tree level results are very close to each other. 

Up to now, we used the common value of the dynamical quark mass M q = 350 MeV motivated 
by the instanton phenomenology. We have to mention that with this value the representation of 
a(t) Eq. (|29|) overestimates the value of the vrA^ sigma term obtained with Eq. 1)27(1 and the latter 
can be reproduced by a small variation of M q . Note that, in the xQSM the dynamical quark mass 
is usually treated as a free parameter varying in the range 350 < M q < 450 MeV. So it is important 
to check whether our results are sensitive to the variation of M q . Certainly, due to Eq. (|44|) and the 



parametric smallness of the term <C 1 the strength of the repulsive core is proportional to 

~ 1/Mq. Interestingly, the strength of the intermediate range attraction actually is insensitive to 
variations of quark mass. In Fig.|2]we show the results for M q ~ 410 MeV (dash-dash-dotted curve) 
which reproduces the same value of the sigma term ovjv as Eq. (|27j) with M q = 350 MeV. The dip 
region is unaffected by these variations with a small coherent shift towards shorter distances. The 
same weak quark mass dependence is found in numerical calculations. Summarizing, the soft pion 
limit provides an accurate representation of the central CrTPE force, compares well with direct 
numerical calculations and accounts for all details of the tree level and unitary forces obtained in 
the previous section. 



D. The spin-orbit CrTPE force 

In coordinate space the spin-orbit part of the scalar-isoscalar CrTPE force is related to the 
central potential by Eq. (|17|) . According to the large- N c analysis given in section IV, the LS 
interaction should be much smaller than the central part because it is suppressed relative to Vc(t) 
by ~ V-^c- O ur numerical results for the LS potential are shown in Fig. |SJ Indeed, the r-space 
behavior of Vz,s(r) with the tree level tttt scattering amplitude (dashed curve in Fig. [HJ) is formally 
similar to the central potential (see the insert of Fig.|2J). But contrary to the Vc(r), we find a sizable 
but much weaker repulsive core ~ 80 MeV (insert of Fig. |SJ) and a very small attraction ~ 0.5 MeV 
at intermediate distances r > 1.2 fm. The effect of the unitarization of the tttt scattering amplitude 
(solid curve in Fig. |5J) plays here the same role as for the central potential and makes the LS force 
softer reducing the strength of the potential by the factor two. Note, that in OBE models, the LS 
interaction induced by effective a exchange generates an attractive spin-orbit force. In our case the 
effect is opposite and leads to a repulsive LS potential in coordinate space. In Ref. Q the spin- 
orbit potential was not calculated (though it could be easily done) and hence direct comparison is 
not possible. We only note, that the results should be similar because of Eq. (|T7j) and because of 
the similar behavior of Vc(r) in both approaches. 

At this level the comparison with other models is also useful. For example, in the standard 



Skyrme model 



321 ] . which includes the non- linear a model Eq. and a stabilizing fourth-order 



term, the "wrong sing" of the isospin independent spin-orbit potential is a long standing prob- 
lem |63j, |64 1 . Several attempts have already been done in this framework to generate the needed 
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FIG. 5: The spin-orbit part of the CrTPE potential in coordinate space. Notations for the curves are the 
same as in Fig. |5J The insert shows the entire structure of the LS potential. 



attractive spin-orbit force, extending the standard Skyrme Lagrangian by including sixth-order 



651 ] spin- 



derivative terms. But the final results differ, varying from repulsive 64] to attractive 
orbit potential. Our results support the existence of a repulsive LS interaction in the isoscalar 
channel. 



E. Analytical formulae for the spin-orbit force 



Here we show an analytical representation of the spin-orbit potential in coordinate space in 
terms of elementary functions. Using Eqs. Q17|) and Q39JI . the LS part Vj^l, which is non- vanishing 
in the chiral limit, reads 

J K 



vS(r) 



64\Z2tt 2 



R S M 2 N 



2304 + 576 £ + 816 £ 2 + 100 £ 3 - 6 £ 4 - 9 £ 5 
2 + iY (A + i+2J2i) (a + |-2 



(47) 
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FIG. 6: The spin-orbit part of the CrTPE potential in the soft pion limit. Notations for the curves are 
the same as in Fig. 0] 

(2) 

Its symmetry breaking counterpart V L $ vanishes in the chiral limit, M v — > 0, and is obtained 
simply from Eqs. (JUJ), (JUiJ) and (glj) 



V$(r) 



J_1 V (2) 

2M 2 r 2 c 



(48) 



4 -3£ 



r 2 \ 3 ) M% 

where for £ = r 2 /i?| < 4 and £ > 4 the V^(r) is given by Eqs. 1)40(1 and ([41 1). respectively. The 



16 + £ 2 2 +e 



Taylor expansion around r — ► up to 0(r 4 ) reads 



r->0 



v/2 i? s M 2 



(49) 



N 



15^vr 2 _^ 

8 i?3 M 2 r + U ' ^ ' 



5^2 



7T~ 



36 



L g 1V± N 

Ml 



S M 2 



(50) 



A' 



40 



i? s M 2 
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and defines the short-range behavior of the LS force with the total strength at origin 



Vls{i 



r^o 3 \ 10 12 * s J R„ M% 



+ ^Xhnj2' (51) 



In the asymptotic region, r —* oo, the soft two-pion tail of the LS" potential behaves as follows 
VS(0| = -192^/X^ + lfe- 10 ) (52) 

VSW =^A T §i + 0(0 (53) 



16V2tt 2 2 7 Mf 1 
3 J " s M 2 , r 6 



The result of Eqs. (|47j) and (|48|1 with the quark mass M g = 350 MeV are shown in Fig. EJ The 
soft pion limit well reproduces the generic features of the LS force (solid line) and compares well 
with our numerical results (dot-dashed curve in Fig. [HJ). Contrary to the central force, even smaller 
contribution of the chiral part V L g{r) (dot-dashed curve in Fig. |EJ) with repulsive strength at r = 
of about ~ 0.11 MeV is noted. The term v£2(r) non- vanishing in the chiral limit accurately 
represents the behavior of the LS potential in all ranges of distances and is actually insensitive 
to the variation of quark masses in the intermediate region. The increase of the quark mass to 
M q = 410 MeV leads to a small shift of the attractive dip towards shorter distances (dash-dash- 
dotted curve in Fig.EJ). At short separation scales, Eq. ([51 j). M 2 R 2 S / '12 -C 1 and the effect of quark 
mass is linear ~ M q and opposite to the central force where the strength is driven by ~ 1/M q . 



F. The ratio V LS /V C and large- N c QCD scaling 

According to large- iV c counting rules discussed in section IV and Eq. ([20(1 . one can expect that 
with the actual number of colors N c = 3 the ratio Vls/Vc should scale according to the following 
simple relation 

Vls/Vc ~ 0{l/N 2 c ) = 0(1/9) ~ 0.11 (54) 

Our results support these expectations and Eq. (|54jl for the relative strength of scalar-isoscalar NN 
components, 1Z = Vls , (0)/Vc(0), is in a good agreement with numerical calculations. With a(t) 
defined by Eq. (|27|) and with the tree level tttt correlation function the ratio is: 1Z ~ 1/7.26 ~ 1.14. 
The unitarization of the tttt scattering amplitude, Eq. ()30|) . plays here an important role and 
improves this value to 1Z ~ 1/9.2 ~ 0.11. Considering the latter value as our parameter free 
prediction, it is remarkable that such detailed information about the relative strength of potential 
components can be deduced directly from QCD. Note that the l/N c expansion is equivalent to 
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FIG. 7: The scaling of the central CrTPE force Vc(r) relative to the spin-orbit potential Vis(r). The 
numerical results are shown with the tree level (dashed curve) and unitary (solid curve) nn correlation 
function. The gray area is the soft pion limit with quark masses M q = 350 (upper limit) to 410 MeV (lower 
limit). The dot-dashed curve corresponds the soft pion limit but with R s related to the qa- The dashed 
horizontal line represents the large- N c QCD expectation, 1/72. = 9. 

the short-distance expansion 2j fif>]. This is the reason, why Eq. (|54|) so accurately represents 
the short-distance part of the CrTPE force. At the same time, it is very difficult to argue that 



381 ] . and 



the scaling of scalar-isoscalar nuclear forces, first predicted by Kaplan and Manohar 
the observation which we made here, are not accidental, because it is very difficult to justify an 
phenomenological models at such separation scales where the baryon number is not well defined 
67] and where other hard QCD processes 68] should play some role. We also recall that the xQSM 
itself as provided by the action Eq. (|21j) is valid for the values of the quark momenta up to the 
UV cut-off A = l/(p) — 600 MeV or r > r\ ~ 0.3 fm. In addition, we have to take into account 
the limited accuracy of the method used here to calcu ate cr(t). Considering our findings as an 
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, should be employed here to check our 



empirical fact, we note that, other models, e.g. Ref. 
statements. To further support the scaling relation Eq. (|54j) we show in Fig. [7| the r-dependence 
of the inverse ratio l/lZ{r) = Vc(r)/VLs( r )- Indeed, the ratio Vc(r)/VLs( r ) scales according to 
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Eq. (|51)l forming a plateau up to distances r ~ 0.8 fm and fastly decreases beyond. The plateau 
covers the region where the action Eq. (j21|) is already reliable and the present approach is well 
justified. 

A good agreement with the scaling relation is also found for the soft pion limit. For M q = 
350 MeV and 410 MeV, Eqs. (jSJ and ((ST) result in Tl ~ 0.12 and ~ 0.16, respectively. The 
radial dependencies of the inverse ratios are shown in Fig. [7| by the band, where the upper limit 
corresponds to M q = 350 MeV and lower one to M q = 410 MeV. Additionally, due to the smallness 
of M%B% in Eqs. (gU) and (ED, the following relation holds: K ~ (9/10)(M% /Mfy. Using Eq. pijl 
with R s = 1/Mq it can be written alternatively 

With the empirical value for the axial coupling constant, qa — 1-25, Eq. (|55|) results in 1Z ~ 0.12 
in good agreement with Eq. (|54[). The r-dependence of 1/7?. with the soliton size i? s related to the 
axial coupling constant g^, Eq. (|24j). is shown in Fig. [TJby the dot-dashed curve. 

VII. SUMMARY 

In summary, we have considered the effect of correlated two-pion exchange modes on the central 
and spin-orbit parts of the scalar-isoscalar NN interaction. We have coupled the tttt correlation 
function to the scalar form factor of the nucleon which we calculated in the large-A c limit in 
the framework of xQSM. The CrTPE force obtained here confirms an unconventional behavior of 
scalar-isoscalar tttt correlations in the NN interaction recently found in unitarized 
the central NN potential with, both, the tree level and unitary tttt correlation functions we find 
a strong repulsive core at short ranges and a moderate attraction at intermediate distances. This 
result indicates that strong repulsive interactions can be generated by the pion-pion dynamics and 
scalar-quark densities themselves. The strength of the intermediate range attraction is insensitive 
to the variation of quark masses. The long-range tail of the central CrTPE AW potential is driven 
by the the ttN sigma term and consistent with the effective a meson exchange. In addition, we 
find a sizable and repulsive spin-orbit force which differs from the LS interaction generated by the 
effective a meson in OBE models. 

The large--/V c behavior of the CrTPE potential was considered and consistency with large-A c 
QCD expectations was found. Both, the central and LS forces satisfy large--/V c QCD counting 
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26J. For 



rules. We have shown that the spin-orbit part is 0(1/ N^) in strength relative to the central force 
resulting in the ratio ~ 1/9 suggested by the 1/N C expansion for N c = 3. The latter is in agreement 
with our numerical analysis and with the Kaplan-Manohar large- iV c power counting. 

Analytical representations for the CrTPE forces in the soft pion limit in terms of elementary 
functions were derived and their chiral content was studied. It was shown that the latter consists 
of two terms. The first one is generated by the symmetry breaking part of the mesonic Lagrangian 
and vanishes in the chiral limit. The second one, which is non- vanishing in the chiral limit, does 
not contain any pion mass dependence and appears as a dominant force which drives all essential 
features of the CrTPE - the strong repulsion at short separations and the moderate attraction at 
intermediate distances. We also find that the soft pion limit provides an accurate representation 
of the CrTPE and compares well with direct numerical calculations. 

In a forthcoming work the role played by the CrTPE in the NN interaction will be addressed 
again and it will be shown that the behavior found here and in Ref. [2| is an important ingredient 
to the NN interaction and appears to be clearly seen in peripheral NN phase shifts above the 
inelastic threshold. 

We acknowledge the correspondence with J.R. Pelaez. This work was supported by the Landes- 
forschungschwerpunkt Baden- Wiirttemberg and the Bundesministerium fur Bildung und Forschung 
(06TU201). 
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